Spin accumulation from the non-Abelian Aharonov-Bohm effect 
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Recently, it has been shown that the spin-orbit coupling (SOC) of the Dresselhaus type in [110] 
quantum wells can be mathematically removed by a non-Abelian gauge transformation. In the 
presence of an additional uniform magnetic field, such a non-Abelian gauge flux leads to a spin 
accumulation at the edges of the sample, where the relative sign of the spin accumulation between 
the edges can be tuned by the sign of the Dresselhaus SOC constant. Our prediction can be tested 
by Kerr measurements within the available experimental sensitivities. 
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Spin transport in the presence of spin-orbit coupling 
has attracted great attention in the field of spintronics. 
The theoretical prediction of the spin Hall effect [3, Q 
in both p and n type doped semiconductors has lead to 
the experimental observation of spin accumulation at the 
sample boundaries 3,4]. Spin-orbit coupling is the main 
cause of spin decoherence in solids. However, recently, 
it has been realized that certain types of spin-orbit cou- 
pling, including models with equal Rashba and Dressel- 
haus SOCs and the model with Dresselhaus SOC in [110] 
quantum wells, can be mathematically removed by a non- 
Abelian gauge transformation In these systems, the 
spin life time is rendered infinite at a magic wave vector, 
giving rise to the Persistent Spin Helix [5, 6]. 

In this letter, we show that the equivalence of the Dres- 
selhaus [110] SOC with a pure non-Abelian gauge flux has 
another intriguing and observable physical consequence. 
In the classic understanding of the integer quantum Hall 
effect, Laughlin and Halperin considered a Gedanken ex- 
periment 0, Q in which one adiabatically inserts a pure 
gauge flux through a cylinder, where the two-dimensional 
electron gas (2DEG) exists on the cylindrical surface, see 
Fig. [U To our knowledge, such a Gedanken proposal 
has never been realized experimentally. We shall show 
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FIG. 1: (Color online) Schematic picture of Laughlin- 
Halperin gauge argument in the Dresselhaus [110] model. The 
Dresselhaus SOC is equivalent to two pure gauge fluxes, act- 
ing oppositely on each spin orientation. The magnetic field 
points out of the cylindrical surface where the 2DEG lives. 
Double arrows indicate the direction of the Lorentz force, 
which leads to the spin accumulation on the edges of the cylin- 
der. 



that the Dresselhaus [110] SOC provides an exact phys- 
ical realization of this Gedanken experiment, where the 
gauge flux is opposite for the two different spin orienta- 
tions. The Dresselhaus SOC can be tuned continuously 
by changing the thickness of the [110] quantum wells 
(QWs). By the Faraday's law of induction, applied sep- 
arately to both spin orientations, the adiabatic turning 
of the gauge flux, or equivalently, the Dresselhaus SOC 
3, leads to a non-Abelian electric field along x direction 
It, £ x = 2md t (3(t)o- z /h 2 . The non-Abelian electric field 
drives electrons with different out-of-plane spin compo- 
nents to opposite y (Hall) directions, resulting in the spin 
accumulation at the sample boundaries. Furthermore the 
relative sign of the spin accumulation between the edges 
can be tuned by the sign of the SOC constant, and the 
spin accumulation has a periodic dependence on the SOC 
constant, with the equivalent Aharonov-Bohm periodic- 
ity. This effect can be detected by Kerr measurements 
within the available experimental sensitivities. 

Our starting point is a 2DEG system with Dresselhaus 
[110] SOC interactions being turned on adiabatically at 
some initial time, and a small Rashba term which we 
treat perturbatively. The electrons are confined in the 
x-y plane subjecting to a magnetic field B = Bz, B > 0. 
There are two intrinsic length scales in our systems, 
namely the magnetic length 4 = ^Jh/eB and the char- 
acteristic length of spin-orbit interaction l so = H 2 /mj3. 
For a fixed magnetic field the ratio r = lb /I so plays the 
role of dimensionless SOC constant. We take the periodic 
boundary condition in x direction and the open bound- 
ary condition in the y direction. In the Landau gauge, 
A = —Byx, the single-particle Hamiltonian is given by 
H{r) = H (r)+H R with 



H = 



(P 



+ eA) 2 2(3 g s y B B 

— -(p x + eA x )a z + ^-a z + V, (1) 



where (3 is the Dresselhaus SOC constant, — e, m and 
Hb are respectively the charge, effective mass and Bohr 
magneton of the electron, the lateral confining potential 
V(y) vanishes for y 6 (— L y /2, L y /2), and is infinite oth- 
erwise. In the following, we first treat analytically the 
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FIG. 2: (Color online) v na as a function of y\ = (±L y /2 — 
Vo)/h, the orbital center measured from the boundaries in 
units of h, for positive ('+') and negative (' — ') edges. The 
energy spectrum at r = is shown in black, and those at 
r = 0.2 for up and down spin are shown in red and blue. 



system Hq{t) and then discuss the effect of the Rashba 
coupling perturbatively. 

For the system Hq we observe that good quantum num- 
bers are the eigenvalues p x — hk and a z = <r = ±1, 
therefore Hq is diagonalized automatically and we seek 



solutions of the form ip(x, y, r) 



,,tkx 



4>nk(y,r) where 



(j> n k is a two-component spinor, obeying the one dimen- 
sional (Id) Schrodinger equation 
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In the above yo a = yo — 2lhar is the orbital center 
for the spin component a and yo = l\k is that with- 
out SOC. We notice that the (3 term in Eq. (1) cor- 
responds to a pure gauge flux since it can be elimi- 
nated by a local non-Abelian gauge transformation U = 
exp (—i2m/3xa z j h ) with the magnitude of the flux given 
by <&<t = 4:in(3L x /h 2 cx r for f3 > 0. Therefore, when in- 
creasing the flux $ CT or r adiabatically, the orbital centers 
of up and down spin components move to the negative 
and positive Hall directions respectively, leading to ac- 
cumulations near the opposite edges. The eigenvalues of 
this system are given by E nka = e nka - hw c (2r 2 - g*a/A) 
where lj c — eB/m and g* = g s m/m e . If we denote the 
Fermi energy of our system as Ep = hw c VF which de- 
pends on the electron density n e and the magnetic field 
the effective Fermi levels for spin-cr are 



° ~ e(i/fr+l/2) ' 

then given by vp a = vp + 2r 2 — 1/2 — g*a/4. We assume 
that the width L y between the two boundaries is large 
enough so that the two edges are well-separated, then in 
the bulk region where \yo a ± L y /2\ 3> h, the edge effect 
can be neglected and £ n ka = huj c (n + 1/2) with integers 
n. While in the edge region where \yo<j ±L y /2\ ~ 4, the 
confining potential must be taken into account by setting 
4>nka{y = TL y /2,r) = 0. In such a case if we still keep 
the form e rikcr = tkjL> c [i/ na (k) + 1/2], however, v na is not 
necessary integers any more, the general solution [l(| to 
Eq. (2) may be written as 



<t>nka{y,r) =e 
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where F[a; b; z] is the confluent hypergeometric function 
of the first kind. The constants C and D are fixed by 
requiring <f> n ka{y - > ±oo) = when considering the 
negative and positive edges separately, which leads to 



r(i+iW2) 



as well as the normaliza- 



£=±2tan(^) 

tion of the wavefunction. The discrete eigenvalue spec- 
trum at each given r are then determined by finding the 
zeros of 4>nka{.y = TL y /2,r). The energy spectrum at 
r = and 0.2 are shown in Fig. [2l 

To discuss the spin polarization, edge charge and 
spin current as well as the Hall conductance, we'll fo- 
cus below on their general dependence on the magnetic 
field and the adiabatical change of SOC constant, com- 
pared in particular with the well-studied case of r — 
0. Throughout this paper we denote the correspond- 
ing quantities in the consideration of the positive (nega- 
tive) edge by a subscript '+(—)'; under such a convention 
the out-of-plane spin polarization near each edge gives 
(a z {r))= Jdyj:Z + \^ nk+ (y,r)\ 2 - YZk~ \^nk-(y,r)\ 2 , 
where y^> = (j/=FL y /2)/4 G [0, Too) is the position vari- 
able measured from the positive (negative) edge in units 
of magnetic length. We emphasize that due to the ex- 
istence of the edges, when vp a exceeds the edge energy 
for some LL n, plays the role of the Fermi 

level for this LL instead, and the Fermi level difference 
between the spin components is no longer g*/2, (which 
is small in general) but increases greatly. This feature 
ensures the large enhancement of the spin polarization 
near the edges relative to the initial (equilibrium) case. 

For r — we reproduce the well-known results for the 
usual 2DEG in which the electrons are magnetized by 
the Zeeman interaction with only very small negative 
spin polarization accumulated near both edges. In this 
case the spin polarization oscillates with 1/B with almost 
constant amplitude. This oscillation pattern is nothing 
but the Shubnikov-de Hass (SdH) oscillation with the pe- 
riod AsdH(l/-B) = lA'&o^e) where $ = h/e is the flux 
quantum. For n e = 5 x 10 11 cm~ 2 [llj], our analytic 
calculation gives a period of 0.048 (T _1 ) which is in pre- 
cise agreement with our numerical data. Interestingly, 
when the SOC constant, positive for example, is turned 
on adiabatically, although the SdH oscillation remains 
with the same period, its amplitude is no longer con- 
stant but changes in different way with 1/B for opposite 
edges. Near the negative edge the amplitude increases 
from negative to positive with the decreasing of the mag- 
netic field, while near the positive edge the amplitude not 
only remains negative but also decreases further. More- 
over, as we further increase r, more up (down) spins are 
pumped to the negative (positive) edge due to the larger 
non-Abelian gauge flux. As an example, the spin polar- 
ization for r = 0.2 and 0.4 are shown in color in Fig. 
[3Ja), among which the results in red correspond to the 
experimental parameters for AlGaAs [110] QWs ll|. To 
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FIG. 3: (Color online) (a) Spin polarization as a function 
of 1/B. The results for r = 0,0.2,0.4 are shown in black, 
red and blue. Insets: Perturbed spin polarizations for r = 
g* = 0.2 at 8 = 0.2, 0.4, 0.8 are shown in green, pink and 
cyan, (b) Spin polarization as a function of r. The results for 
vf = 1.0,2.0,3.0 are shown in black, red and blue. In both 
figures, right panel: negative edge; left panel: positive edge, 
n e — 5 x 10 n cm~ 2 and the solid, dash-dot lines indicate the 
results for g* = 0.2, 0.4. 



have an overview of the dependence of the spin polar- 
ization on the SOC constant, we also plot in Fig. [3Jb) 
the spin polarization as a function of r at the Fermi en- 
ergies vp = 1.0,2.0,3.0 respectively, the result of which 
shows clearly that the larger the SOC constant is, the 
more positive (negative) the spin polarization near the 
negative (positive) edge is. Furthermore, there is a simi- 
lar oscillation pattern with increasing r, where the turn- 
ing point r„ occurs whenever the gauge flux changes to 
the values so that the Fermi energy for spin a crosses a 
bulk LL vp a = n, and the period of which is given by 
Ar 2 (<f> cr )| r<T = 1/2. This periodicity corresponds exactly 
to the Aharonov-Bohm period for both down and up spin 
components of the gauge flux. For negative SOC constant 
the same results are obtained with only the sign-exchange 
of the spin polarization between the boundaries. The 
above results present an interesting scenario for the spin 
accumulation: electrons with opposite out-of-plane spin 
components accumulate near the opposite edges, the rel- 
ative sign of which can be tuned by the sign of the SOC 
constant; by decreasing the magnetic field or increasing 
the SOC constant and/or the electron density, the mag- 
nitude of the spin polarization can be enlarged greatly 
compared to the equilibrium case. 

This phenomenon is a direct result of the intrinsic elec- 
tric SU(2) gauge field in our system. Mathematically it 
is clearly seen through the expression of j/oo- that for fixed 
lb the adiabatical increasing of r (positive) moves the or- 



bital centers of up spin to negative y direction, and in the 
meanwhile carries those of down spin to positive y direc- 
tion. The same conclusion is reached through similar 
analysis when decreasing the magnetic field for fixed r. 
To have an idea of the magnitude of the spin accumula- 
tion in real materials, we recall that the range of the SOC 
constants from weak to strong is [5j 10~ 13 — > 10~ n eV-m 
corresponding approximately tor £ (0, 0.5) which is cov- 
ered by our results shown in Fig. [3J 

Let's now turn our attention to the charge and spin 
current as well as their Hall conductance. The charge 
current operator for each spin component is given by 
j xa = —ev xa - where v xa — d Q° a is the velocity opera- 
tor. The total charge current is the sum of that from up 
and down spin components J£(r) = eui c h J ^ES©"'" 
yo + 2ra)\cj) n ka\ 2 , and the charge Hall conductance (CHC) 



follows G% y = G c xy:+ 



G%y,-- To compare we also 



consider the spin current polarized in z direction flow- 
ing along the edges by taking the traditional definition 
3x = j{ v x,crz} = — {fr/2e)jx a zi where we see that in our 
system this spin current is physically just the ovcarrying 
charge current. Therefore the definition of spin current is 
physically meaningful here in the sense that it satisfies a 
well-behaved continuity equation since the charge current 
does, so there is no controversy on the conservation issue 
which usually can't be avoided in general SOC systems. 
Different from that of charge current, the total spin cur- 
rent is given by the difference of the currents from up and 
down spin components 1%. (r) = J* , (r) — 1*_ (r) and so as 
the spin Hall conductance (SHC) G s xy — G xy + — G s xy _. 
Our numerical results for CHC and SHC at the param- 
eters r = g* = 0.2 close to those for AlGaAs [110] QWs 
1 1( | are shown in Fig. [4] it is seen that for CHC there are 



wider steps when it is quantized in even multiples of the 
quantum of conductance, while the steps are much nar- 
rower when it is quantized in odd multiples of the quan- 
tum of conductance. For SHC they are non-vanishing 
only when the CHC takes odd integers and the corre- 
sponding quantized number is unity. 

The above results can be easily understood by the stan- 
dard argument [lOJ that the charge current flows only 
if there exists a chemical potential difference between 
the edges, and whenever the current flows CHC is quan- 
tized in units of e 2 /h with the quantized number be- 
ing the number of edge states occupied. In our system, 
for each spin component we have G xy a = n + 1 with 

vfo S (Vn- e > u n+i an d the total CHC takes even or 
odd integers depending on the relative fillings of the LLs 
by up and down spin. The quantization of SHC is a 
direct result of the quantization of the CHC of a spin- 
polarized LL in the presence of the Zeeman splitting at 
odd integer filling factor. It is obvious from the rela- 
tion between the charge and spin current operators that 
the SHC for each spin component is quantized in units 
of — e/47r with the quantized number equals to that of 
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FIG. 4: (Color online) Quantized charge and spin Hall con- 
ductance as a function of 1/B. The CHC and SHC are shown 
respectively as black and red lines. The parameters are taken 
as r = 0.2, g* = 0.2, n e = 5 x 10 n cm" 2 . 

the CHC, hence it is nonvanishing only when the total 
CHC takes odd integers otherwise different spin compo- 
nents will cancel each other giving zero. The width of 
the odd steps is A(l/B) = g*AsdH~/2 which equals ap- 
proximately to 0.0048 (T _1 ) for the parameters taken in 
Fig. [4] It is worth noting that in a previous work by Bao 
et al. [l2T |. similar results for charge and (resonant) spin 
Hall conductance are first obtained in a general model 
with both Rashba and Dresselhaus SOC, where the SHC 
is not quantized and the spin accumulation decays due 
to the diffusive property in these systems. 

Finally we discuss the Rashba-type perturbation term 
Hr = 5/3{p y <7 x — (hk — eBy)<Ty]/h, where the Rashba SOC 
constant has been written as a = 8(3 with 8 being a small 
perturbation number. This Hamiltonian couples the LLs 
of different spin components and the wavefunction per- 
turbed to the first order of Hr gives 

t , , \ ^ t^n' k&'.nka , / ,\ 

Qnka = <Pnka + ^ ^ -pn'ka , (4) 

n' 

where E nka - E n >kw = hu c [v n<y - v n <-a + g*(a -a)/A] 
and Wn'ka-nka = J dy(j)* n , Kff(y)H R,wa(y \ dy)4>nka{y) with 
a = —a. We have studied systematically the effect of the 
perturbations to the spin polarization and some of the 
results are shown as insets in Fig. [3ja). It is concluded 
that the perturbations in general make the spin polar- 
ization more positive at both edges so that (cr z (r)) is in- 
creased near negative edge while decreased near positive 
edge, and the larger 8 and/or r is the more the increases 
(decreases) are. But in general the magnitude of these 
deviations is very small and will not qualitatively change 
our picture of the spin accumulation. 

We propose to use Kerr measurements to detect our 



predictions on spin accumulation. It is calibrated in the 
experiment [3j that for the film-like 3d sample, the spin 
density of 20 Bohr magnetons per /zm -3 is signaled by 
1 fiv&d of Kerr rotation angle. Bearing in mind that the 
film thickness used in the experiment I = 0.6 /zm is 2 to 
3 orders smaller than the other two dimensions, hence 
a similar correspondence can be estimated for quasi-2d 
systems that the spin density of 10 Bohr magnetons per 
/im~ 2 lead to 1 /irad of Kerr rotation angle. For the 
parameters of AlGaAs [110] QWs [ll| used in our cal- 
culation, the spin density n s = (tJ z )/l% is estimated to 
be 10 4 Bohr magnetons per /zm -2 at B — 5 T, which 
is several orders above the state-of-art detection limit in 
the experiment [H, and our predicted spin accumulation 
can be easily detected within the available experimental 
sensitivities. 
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